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Abstract 

The isoscalar El and E2 resonances in the proton-rich nuclei ^^'^^Ni and the {/7/2<^ 
3~} multiplet in ^^Ni have been calculated taking into account the single-particle con- 
tinuum exactly. The analogous calculations for the mirror nuclei ^'^Ca and ^^Sc are 
presented. The models used are the continuum RPA for ^'^Ni, ^^Ca and the Odd RPA 
for ^^Ni, ^^Sc, the latter has been developed recently and describes both single-particle 
and collective excitations of an odd nucleus on a common basis. In all four nuclei we 
obtained a distinct splitting of the isoscalar El resonance into Ihuj and 3huj peaks at 
about 11 MeV and 30 MeV, respectively. The main part of the isoscalar El energy- 
weighted sum rule (EWSR) is exhausted by the 3hu! resonances. The IHu! resonances 
exhaust about 35 % of this EWSR in 48,49]vgi ^^^^^^ 22 % in ^^Ca and ^^Sc. Ah 
seven {/7/2 "X) 3~} multiplet members in ^^Ni are calculated to be in the (6-8) MeV 
energy region and have noticeable escape widths. 



PACS numbers: 21.60.Jz, 24.30.Cz, 27.40. +z 
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The main feature of the proton-rich unstable nuclei ^^Ni and ^^Ni, which have been 
discovered recently [||, 0, is that they have a proton separation energy near zero. It seems 
to be of great interest to calculate excitations of these nuclei, especially in the framework of 
the same calculational scheme, which would take into account the single-particle continuum 
reliably, because its role is obviously very important in this case, though the corresponding 
effects may be not so strong as in the neutron-rich nuclei with the neutron separation energy 
near zero. The comparison of the theoretical predictions with future measurements of the 
excitations of these nuclei will make it possible to obtain an important information about 
the nucleon-nucleon effective interaction in drip-line nuclei and about features of such nuclei, 
which are also of astrophysical interest. 

For the doubly-magic ^^Ni, one of the proper microscopic methods is the theory 
of finite Fermi-systems (TFFS) |^ with exact accounting for the single-particle continuum 
(continuum TFFS). It corresponds to the continuum random phase approximation (RPA) 
making use of the phenomenological Landau-Migdal particle-hole (p-h) interaction and a 
single-particle Woods-Saxon scheme. The parameters describing this interaction and the 
single-particle scheme are known 0, ^, |^ , at least for stable nuclei, and it is of great interest 
to check them for unstable ones. 

For doubly magic nuclei plus or minus one nucleon, a theory has been developed 
recently 0, which was termed the Odd RPA (ORPA). It corresponds to the continuum RPA 
for even-even nuclei and, in addition, consistently accounts for properties of odd nuclei. The 
ORPA allows a consistent description of both the single-particle and collective part of the 
excitation spectrum of an odd nucleus without pairing, including giant resonances in the 
continuum and splitting of particle (hole) ® phonon multiplets, on a common basis. The 
description of these multiplets taking the single-particle continuum into account is especially 
important and interesting for nuclei such as ^^Ni with the nucleon separation energy near 
zero and relatively large phonon energies. Other questions concerning the excitations of this 
nucleus exist. For example, it was shown in that the isoscalar (IS) El strength in ^''O 
below 12.5 MeV is completely determined by the odd nucleon contribution. The question 
about a similar effect in "^^Ni arises. 

Thus, we have a reasonable and very interesting possibility to calculate excitations 
of both nuclei within the same calculational scheme based on the continuum TFFS as well 
as to investigate primarily the role of the small proton separation energy in the proton-rich 
48,49]\jj nuclei and of the single-particle continuum. This is the aim of the present article. 
We will present the results of the calculations of the IS El and E2 giant resonances in ^^Ni 
and ^^Ni. For ^^Ni the characteristics of the {/7/2 ® 3~} multiplet will be presented for all 
7 resonances. For comparison and reliability of the results, the resonances in the mirror 
nuclei ^^Ca and "^^Sc and the multiplet in ^^Sc are also calculated. In the calculations we 
use the TFFS and Woods-Saxon potential parameters which are mainly known. However, 
in order to have reliable results for the exotic nuclei, we fit some of the parameters of the 
Landau-Migdal interaction and Woods-Saxon potential using Hartree-Fock single-particle 
levels obtained with the known Skyrme interaction and experimental data for ^^Ca nucleus 
(see below). 

In the TFFS approach the effective interaction and the single-particle scheme 
are described by two sets of phenomenological parameters and are not connected with each 
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other. In our calculations the single-particle potential was taken in the form 



U,{r) = Ulfwsir) + U^^^^^^il ■ a) + Uc{r) , (1) 

where fws{r) = 1/[1 + exp((r — R)/a)], index q = n,p denotes the sort of the nucleon. The 
Coulomb potential Uc for protons is that of a uniformly charged sphere. The parameters 
Ug°, n, R, and a used here are well known P: 

Vpp = Vnn = l9-7 MeV , V^^ = V^p = 87.0 MeV , 
K = 0.263 (1 + 2{N - Z) I A) fm^ , 
i? = 1.24A^/3fm, a = 0.63fm. 

The choice of the parameters in Eq. (|l|) is a more complicated question for the 
48,49]\j-j nijicigi xhe point is that the value of cannot be determined from the standard set 
of parameters @, which corresponds to = in Eq. (|l|), because in this case it leads to 
a positive energy of the last occupied proton I/7/2 level in "^^Ni. To overcome this difficulty 
we used the following method. We have calculated the energies of single-particle proton and 
neutron I/7/2 levels in ^^Ni in the framework of the self-consistent Hartree-Fock method with 
the Skyrme SIII interaction of Ref. which is often used in the calculations of properties 
of both stable and unstable nuclei. Then the values of the potential depths and f/° in 
Eq. (|1]) were fitted so that to reproduce the Hartree-Fock values of £^(1/7/2) = —0.86 MeV 
and e„(l/7/2) = -16.98 MeV. From this input we obtained [/° = -52.03 MeV, f/° = -60.98 
MeV for ^*Ni. The same procedure was performed for the ^^Ca nucleus. In this case the 
values = -58.32 MeV and t/° = -49.09 MeV yield the Hartree-Fock values of the single- 
particle energies £p(l/7/2) = —10.16 MeV and £:„(l/7/2) = —9.99 MeV. The corresponding 
experimental values for ''^Ca are —9.63 MeV and —9.95 MeV, i. e. the agreement with 
experiment is fairly good. 

The Landau-Migdal p-h interaction used in our calculation is defined by the fol- 
lowing ansatz 

^(n, r2) = Co (/(ri) + /'(ri)(Ti ■ T2) + g{(T, ■ a,) + g\Ti ■ T2){(t, ■ a^)) 6{r, - r,) , (2) 
where 

/(r) = fe. + {hn - /e.)p(r)/p(0) , f {v) = + {f^ - /L)p(r)/p(0) . (3) 
The total nucleon density p(r) entering Eqs. (|^) is defined in our approach by the formula 

p{y) = ^nx\ipx{r,a)W (4) 

X,cr 

where n\ are the occupation numbers, ip\ are the single-particle wave functions calculated 
within the above- described Woods-Saxon scheme. 

The parameters entering Eq. are known [^, |^, as a rule, except for one of the 
six parameters, namely the parameter f^x, which is fitted very often to the experimental 
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energies of low- lying phonons, see, for example, Refs. 0, ^. In our case this is especially 
necessary because, strictly speaking, there is no guarantee that the Landau-Migdal inter- 
action is suitable for the drip-line nuclei. Since the excitations of ^^Ni are unknown, it is 
reasonable to use known excitations of the mirror nucleus ^^Ca for adjusting f^x- So we have 
fitted the parameter f^^ to reproduce (within the continuum TFFS) the energy of the 3]~ 
level in ^^Ca at 4.51 MeV. The result was fe^ = —0.94. Thus, in all our calculations the 
following set of the Landau-Migdal interaction parameters has been used (the rest of them 
being taken from Refs. [§ |, |, 0) 

/e, = -0.94, /,„ = -0.002, /L = 2.30, = 0.76 , 

g = 0.05 , g' = 0.96 , Co = 300 MeV ■ fm^ . 

With these parameters we obtained a reasonable agreement with experiment for 2^ 
level (^th(2+) = 4.15 MeV, E^^pC^t) = 3.83 MeV) and very good agreement for the B{E3) ] 
= B{E3, g.s. 37) value {Bth{E3) T = 6500 e'^fm^ and B^xp{E3) ] = 6700 e^fm^) in ^^Ca. 
The latter is important for the realistic description of the {/7/2 ® 3~} multiplet. 

The calculations of the IS El and E2 resonances in ^^Ni and ^'^Ca were performed 
within continuum TFFS. For the IS El resonance the "forced consistency" and the spurious 
state suppression procedure was applied. This procedure and other details of the calcula- 
tional scheme have been described in Refs. |0, |1^. The same resonances and the {/7/2 ® 3~} 
multiplets in ^^Ni and ^^Sc have been calculated within the ORPA model developed in [0. In 
the ORPA calculations, only one (entrance) channel was incorporated and same interaction 
parameters and the single-particle schemes were used as in our TFFS ones. All the calcula- 
tions were performed using the coordinate representation technique in the framework of the 



Green function method |TT[ which enables one to take into account the single-particle con- 
tinuum exactly. The smearing parameters used were A = 300 keV for the giant resonances 
and A = 0.2 keV for the multiplets. 

One-body transition operators (i. e. the external field operators in terms of 
TFFS) were taken in the following form: 

V^oiEi) ^ (^3 _ R2^r)Y,,{r) (5) 

for isoscalar El excitations, 



for isoscalar E2 excitations, and 



yOiES) ^ e^Y.^if ) (7) 



for electromagnetic E3 excitations, where Cq are the nucleon effective charges in the center- 
of-mass reference frame (see, e. g., Ref. [pi]]): 



6p 



with L = 3. Parameter Rjs in Eq. (^) is determined within the "forced consistency" scheme 
(see 0). In the present calculations the numerical values of Rjs for all nuclei under consid- 
eration are very close (within 1 %) to those of (| < >)^/^ (in this context, see Ref. |T^). 
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The results for the IS El and E2 resonances are presented in Figs. Q and ^ in the 
form of the normahzed energy-weighted strength functions Sew{E): 

Sew{E) = ^ES{E) (9) 

where S{E) is the RPA (ORPA) strength function, mf^ is the corresponding energy-weighted 
sum rule (EWSR), which is defined as 

^ 3^ / 25 ^2 2\ ^^Q) 

^ 8m7r V 3 y ^ ' 

for the IS El operator (|^), and 

mf^ = ^ < > (11) 
4m7r 

for the IS E2 operator (^. Note that these formulas differ from the often used definitions of 
the IS EL EWSR (see Refs. [|1^, |^) by the factor (2L + 1) arising due to the summation over 



the quantum numbers of the excited states and of an external field in the strength functions. 

The functions Sew{E) for IS El excitations of all the nuclei considered are shown 
in Fig. m. The IS El strength is distributed around two main maxima at about 11 MeV and 
30 MeV, which correspond to the Ifuv and ?)t\jjj IS El resonances. The calculated position 
of the ?)kuj resonances is in agreement with the experimental formula E^^ ~ 114^4^^/'^ = 31.4 
MeV for A = 48, which was obtained in Ref. |]14| for the upper components of the IS El 



resonances in Zr, Sn, and Pb. To estimate the contribution of these resonances into 



the IS El EWSR (p!0| ) the following integrals were calculated 

p{Emin', Emax) = J SEw{.E)dE. (12) 
E 

The contribution of the Ihw resonance {p{lhw)) was defined as p(0; 21) in ^^Ni, ^^Ni (here- 
inafter the values of Emin and E^ax are given in MeV), and as p(0; 19) in ^'^Ca, ^^Sc, where 
the values of 21 MeV and 19 MeV are the energies of minima separating Ihu and Shw 
resonances in these nuclei. For the Sfkj resonance the contribution {p{3hw)) was defined as 
p{21; 40) in ^^Ni, ^^Ni, and as p{19; 40) in ^^Ca, ^^Sc. The results, which are listed in Table |], 
show that the main part of the IS El EWSR is exhausted by the 3huj resonances in all nuclei 
under consideration. The Ihw resonance exhausts about 35 % of the IS El EWSR in '^^'^^Ni 
and about 22 % in ^^Ca and ^^Sc. 

Further, comparing the results for the doubly-magic and neighbouring odd nuclei in 
Fig. |l|, we see that the role of the odd nucleon is rather small in our case. The difference from 
the IS El resonances in ^^O and ^^O (see Ref. |Q) is that in ^^Ni and ^^Sc the single-particle 
resonances are embedded into the region of the core p-h excitations and are only slightly 
distinguishable against their background . 

The IS E2 resonances in ^^Ni, ^^Ni, ^'^Ca, and ^^Sc are shown in Fig. Their 
characteristics are given in Table 0. We obtained a reasonable depletion of the IS E2 
EWSR within the 10-30 MeV interval taken. The positions of the maxima of the calcu- 
lated strength functions are in a reasonable agreement with an empirical formula for the 
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Figure 1: The normalized energy- weighted strength functions Sew{E) for the isoscalar (IS) 
El excitations in ^^Ni (a) and ^'^Sc (b) calculated in the ORPA (solid lines), and the same 
functions for ^^Ni (a) and ^^Ca (b) calculated in the continuum TFFS (dotted hnes). The 
smearing parameter is A = 300 keV. 
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Figure 2: Same as Fig. 1, but for the isoscalar E2 resonance. 
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Table 1: The percentages of the IS El EWSR for the Ihu and 3huj IS El resonances (see 
text for details). 









49Ni 


4«Ca 


^^Sc 




% EWSR 


35 


35 


21 


22 




% EWSR 


59 


60 


73 


73 



Table 2: Energy of the maximum E^, the width at half maximum Fq, and the percentage 
of the IS E2 EWSR in the energy interval 10-30 MeV for IS GQR. 





48Ni 


49Ni 


48Ca 


49Sc 


Eo, MeV 


17.8 


17.6 


16.6 


16.4 


To, MeV 


1.5 


1.8 


1.1 


1.3 


p(IS GQR), % EWSR 


94 


94 


91 


92 



mean energy of the IS giant quadrupole resonance (GQR) Eisgqr ~ 63 A^^^'^ MeV = 17.3 
MeV for A = 48 and with the mean energy E = 16.9 MeV obtained for '^^Ca in [|I^ in the 
lplh+lplh(8)phonon+GSCp/ionon approximation. The present results for IS GQR in ^^Ca 
are also in good agreement with those obtained in Refs. |]16|, within the self-consistent 



continuum RPA (in Ref. |jT6| the maximum of S{E) was obtained at 16.6 MeV, and 87 % of 
the IS E2 EWSR was found within the 14-19 MeV interval). However, in our calculation the 
peak energy of IS GQR in ^^Ni is about 1.8 MeV shifted upwards from the value of 16.0 MeV 
obtained in Ref. . This considerable discrepancy seems to be explained by the difference 
between the single-particle schemes and between the effective forces used in the calculations 
being compared. The case of the ^^Ni nucleus appears to be more sensitive to this difference 
due to the small proton separation energy. 

From comparison of the results presented in Figs. |I| and H one can see that the 
proximity of the proton separation energy to zero in ^^'^^Ni is especially manifested in the 
case of Ihw IS El resonance giving rise to a noticeable increase of its escape width F^ and to 
the larger depletion of the IS EWSR (about 35 %) as compared with ''^Ca and "^^Sc (about 
22 %). The broading of the IS E2 resonances in the unstable nuclei can be also seen from 
the results presented in Table 0. 

Note, however, that our calculations cannot pretend to describe total width of the 
resonances because only two main damping mechanisms (Landau and escape ones) were 
taken into account while the third one, i. e. the spreading-width formation mechanism, was 
not included in the calculations (for such kind of calculations in ^^Ca see Ref. |ll5[| ). 

At last but not the least, in Fig. | we show the strength functions for the {/7/2®3~} 
multiplets in ^^Ni and ^^Sc together with the 3]~ levels in '^^Ni and "^^Ca, respectively. Here 
the neutron I/7/2 single-particle level corresponds to the ground state (g.s.) of the '^^Ni 
nucleus, the proton I/7/2 level corresponds to the g.s. of the ^^Sc. As compared with the 
giant resonances (Figs. 0) we see a very natural and noticeable role of the odd neutron in 
"^^Ni and proton in ^^Sc. All seven multiplet members are around the 3j~ level of the core at 
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Figure 3: The electromagnetic E3 strength functions S{E) for ^^Ni (a) and ^'^Sc (b) calcu- 
lated in the ORPA (solid lines), and the same functions for ^^Ni (a) and "^^Ca (b) calculated 
in the continuum TFFS (dotted lines). The smearing parameter is A = 0.2 keV. 
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Table 3: Characteristics of the septuplet members {/7/2 ® in '^^Ni and ^^Sc. The 

calculated energy shifts 6E from the excitation energies of 3^^ states in "^^Ni {E{3i) = 6.802 
MeV) and ^'^Ca {E{3i) = 4.518 MeV) are shown. The escape widths F^^ for the septuplet 
members in ^^Ni are given. The values of B{E3, j'^ g.s.) are shown in Weisskopf units 
(1 W. u. = 142.6 e^fm'^). 



f 




49Ni 




49 


Sc 




5E 


rT 


B{E3) 


5E 


B{E3) 




(keV) 


(keV) 


(W. u.) 


(keV) 


(W. u.) 


1 + 

2 


279 


33 


11.1 


228 


3.4 


3 + 


-244 


49 


13.9 


-14 


6.7 


!+ 


361 


39 


11.4 


142 


5.1 


f+ 


-455 


52 


14.7 


-197 


10.1 


1+ 


33 


77 


16.2 


5 


6.9 


i!+ 

2 


-401 


46 


14.3 


-235 


10.3 


13 + 
2 


1052 


28 


6.7 


140 


1.0 



6.8 MeV in the case of '^^Ni and at 4.5 MeV in the case of "^^Sc. Since the 3f level in ^^Ni 
lies in the continuum (and consequently is a resonance), each of the multiplet members in 
^^Ni has a noticeable escape width Our approach enables one to calculate the escape 
width immediately because the single-particle continuum is included completely. The value 
of was determined as a width at half maximum of the corresponding strength function 
(note that the smearing parameter A = 0.2 keV in Fig. The calculated value for "^^Ni 
is r^(3j~) = 51 keV. For the multiplets members, the values of T\ energy shifts 5E, and 
B{E3, j'^ — > g.s.) are listed in Table |^. It is of great interest to measure the resonances 
in ^^Ni and their "pure" escape widths, which may be an easier task at present than the 
measurement of the 3j~ level in ^^Ni. 

In conclusion, we have calculated some excitations in ^^Ni, "^^Ni, ^^Ca and ^^Sc 
within the same calculational scheme and with the consistent account for the single-particle 
continuum. Despite the fact that more complex configurations have not been accounted 
for in the calculations, our results expected to be realistic enough to be compared with 
experiment. In all four nuclei considered a distinct splitting of the IS El resonance into 
two peaks at about 11 MeV and 30 MeV has been obtained. The influence of the small 
proton separation energy in ^^'^^Ni has been demonstrated, especially on the escape width 
and a large depletion of the IS El EWSR by the Ihw IS El resonance. For the {/7/2 ® 3^} 
multiplet in "^^Ni, the widths of all the multiplet members are formed by the single-particle 
continuum in our approach. Measurements of these excitations, including escape widths of 
the individual states, will give valuable information about many features of drip-line nuclei. 

The authors are very grateful to Prof. V. R. Brown for her advices and careful 
reading of the manuscript. 
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